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ABSTRACT 

We describe observations of aromatic features at 7.7 and 11.3 /xm in AGN of three types including 
PG, 2MASS and SCR objects. The feature has been demonstrated to originate predominantly from 
star formation. Based on the aromatic-derived star forming luminosity, we find that the far-IR emission 
of AGN can be dominated by either star formation or nuclear emission; the average contribution from 
star formation is around 25% at 70 and 160 /im. The star-forming infrared luminosity functions of the 
three types of AGN are fiatter than that of field galaxies, implying nuclear activity and star formation 
tend to be enhanced together. The star-forming luminosity function is also a function of the strength 
of nuclear activity from normal galaxies to the bright quasars, with luminosity functions becoming 
flatter for more intense nuclear activity. Different types of AGN show different distributions in the 
level of star formation activity, with 2MASS > PG > SCR star formation rates. 
Subject headings: infrared: galaxies - galaxies: active - galaxies: starburst 



1. INTRODUCTION 

The interplay between supermassive black holes 
(SMBHs) and star formation is now recognized as a 
critical ingredient in galaxy evolution, as demonstrated 
by the correlations between the blackhole masses and 
the bul ge properties of their host galaxies (M-a re- 
lation) (iKormendv fc Rich stond 11995.: Magor rian et al.l 
Il998t iGebhardt et al. 2000: iFerrarese fc Merr itt 2000^ 
However, because the star formation rate (SFR) is diffi- 
cult to measure around active galactic nuclei (AGN), we 
are unable to answer basic questions about the interre- 
lations between the two processes: in what star-forming 
environments does AGN activity tend to be triggered? 
Does feedback from one process trigger or quench an- 
other? 

Models that involve the galaxy merging process and 
AGN feedback simula t e the M-a relation successfully 
(e.g. iDi Matteo et al.l l2005( l. The theoretical pic- 
ture of the "cosm ic cycle" of galaxy evolution (e.g. 
iHopkins et al.l[2006h connects galaxy mergers, starbursts 
and nuclear accretion. Galaxy mergers induce gas in- 
flow producing starbursts and obscured quasar activity. 
As the quasar feedback starts to heat and expel the cir- 
cumnuclear medium, the nuclear activity becomes visi- 
ble as optically bright quasars. Eventually, the quasar 
activity and starbursts are terminated as the gas and 
dust are more thoroughly expelled. In this scenario, the 
time histories of the star formation and nuclear accre- 
tion through the merging process are two fundamental 
physical properties underlying many ob s ervations (e.g. 
iGranato et"^ 120041 : ISpringel et al1l2005t IHopkins et al.l 
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l2006l ). However, current observations only provide de- 
tailed understanding of star formation in normal galax- 
ies, not in those dominated by luminous AGN. 

While the near- and mid-IR emission of AGN arise 
from hot and warm d ust heated by nuclear emission (e.g . 
Polletta et afl I2OOOI : iShi et all [2005t iHines et all I2OO6I : 



Jiang et al.l the heating mechanism of the cold 



dust responsible for the far-IR e mission still remains am- 
biguous (Sec Ha as et al.l I2003D. As suggest ed by nu- 
merical simulations ([Chakrabarti et al.l 120061 ). the con- 
tribution of the AGN to the far-IR emission may char- 
acterize different evolutionary stages. Insights into the 
far-IR emission mechanism can also constrain the struc- 
ture of the cir cumnuclear mater ial and its evolution with 
redshift (See iBallantvne et al.l [^006). It is also crit- 
ical to understand the energy b udget of many AGN 
revea l ed in deep IR surveys (e .g. lAlonso-Herrero et al.l 
I2OO6I : iDonlev et al.ll2005l l2007| ). which are faint in the 
optical, or even in X-ray bands, and whose main energy 
output resides at infrared wavelengths. Progress on these 
topics requires the ability to identify the contribution of 
star formation to the IR emission. 

Although the commonly used star-formation tracers 
(the total UV, Ha and IR emission) may be contami- 
nated severely by the nuclear emission, there are several 
alternatives to estimate the SFR in AGN, such as the 
extended UV emission, extended mid-IR emission, and 
narrow metal emission lines. The extended UV emission 
can be observed with high-resolution telescopes such as 
HST. However, due to the large brightness contrast be- 
tween type 1 AGN and the host galaxy in the UV, this 
method is limited to type 2 AGN, and even for them 
the scattered nuclear UV emission may be significant 
([Zakamska et a l. 20061 ). Extended mid-IR emission has 
been us ed to estimate the SF R for nearby Seyfert galax- 
ies (e.g. iMaiohno et anil995f ). Due to the hmited angu- 
lar resolution of infrared telescopes, it becomes difficult 
to resolve the AGN from the circumnuclear star forma- 
tion for objects at z >0.05 (0.5"=500pc). Estimating the 
SFR with narrow metal emission lines is difficult because 
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they are contaminated by the AGN narrow emission hne 
region. In addition, this method suffers from other prob- 
lems, for example, the [O II]A3727 flux of PG quasars 
indicates a very low SFR (|HoI [20051) . which IS inconsis- 
tent with the abundant molecular gas in these objects 
and possibly a result of under- estimating the arnount of 
extinction of the emission line (jSchweitzer 

In this study, we employ the mid-infrared aromatic fea- 
tures to quantify the SFR in AGN host galaxies. These 
feat ures are proniinent at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 
/im (iGillett et al.lfl973l) . They are seen in various Galac- 
tic environments including HII regions, diffuse interstel- 
lar clouds, planetary nebulae, reflection nebulae and pho- 
todisassociation regio ns (PDRs) and in ex tragalactic ob- 
jects (for a review, See lTielens et al.ll999f) . The aromatic 
emission in normal star-forming galaxies is similar to 
that in Galactic star- fo rming regions (e.g. iGenzel et al.l 
119981 IClavel et al]|2000l). with a well understood correla- 
tion to the SFR fe.g. iRoussel et al.ll200ll : iDale fc Heloul 
120021 ). The aromatic features in active galaxies have 
much lower e quivalent width (EW) than in star-for ming 
galaxies (e.g. iRoche et al.lll991l: IClavel et al.ll2000[ ). im- 
plying the destruction of the aromatic carriers by the 
harsh nuclear radiation or the inability of the nuclear 
radiation to excite the aromatic features. Evidence for 
excitation of the aromatic features by star formation 
in active galaxies comes from spatially resolved mid-IR 
spectra of nearby examples, where the observed aromatic 
emission is mainly from the disk (e . g. ICutri et al.l 19841 : 
Desert & Dennefcld"l988|; IVoitlll992t llaurent et al]l2000l : 



Le Floc'h et al. 2001). Various infrared diagnostics have 
been developed based on a correlation of aromatic fea- 
ture strength with star forming activity to discriminate 
the power sources (star formation versus nuclear activity) 
for luminous infrared galaxies ( LIRGs; L/j? > 10^^ L(T)) 
(e.g. [Gcnzcl ct al. 1998; Laure nt et all 120001 : iTran et al.l 
I2OOII : IPeeters et al.l [20041 ) . Direct measurements of the 
aromatic features in a small PG quasar sample have 
been carried out bv lSchweitzer et aIT()2006D to study the 
quasar far-IR emission mechanism. 

In this p aper, we present Spitzer Infrared Spectro- 
graph fIRS: [Houck et ani2004D low-resolution spectra for 
a large sample of AGN. § 2 describes the sample, the data 
reduction, the extraction of the features at 7.7 and 11.3 
fim and the determination of the associated uncertain- 
ties. In § 3, we provide evidence for the star-formation 
excitation of the aromatic feature in these objects. In 
§ 4, we estimate the conversion factor from the aromatic 
flux to the total IR flux. § 5 discusses the origin of AGN 
far-IR emission. In § 6, we construct the luminosity func- 
tion of the SFR in AGN host galaxies and discuss its 
implication for AGN activity. § 7 presents our conclu- 
sions. Throughout this paper, we assume Ho=70 km 
Mpc-i, and 

2. DATA AND ANALYSIS 

2.1. Sample 

The sample in this paper is composed of objects 
derived from three parent samples selected by differ- 
ent tech niques: optically-selected Palomar-Green (PG) 
quasars (I Schmidt fc GreenI [19831); the Two -Mic ron All 
Sky Survey (2MASS) quasars (ICutri et al. 2001*): and 
SCR radio galaxies and quasars ([Spinrad et al., il985D . 
PG quasars are selected at B band to have blue U-B 
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Fig. 1. — The redshift distribution of the three subsamples 
in this study (shaded area) compared to the corresponding parent 
samples for the PG, 2MASS and SCR objects. The insert plots 
show the flux distribution of the subsample (shaded area) and the 
corresponding parent sample for the 3CR and 2MASS objects. 



color, a dominant starlike appearance, and broad emis- 
sion lines. 2MASS quasars represent a much redder near- 
IR-to-optical quasar population comp ared to PG quasar s 
but have similar i^s-band luminosity (jSmith et al.ll2003 ). 
Unlike PG quasars, the 2MASS and 3CR samples include 
objects with narrow, intermediate and broad emission 
lines. 

Besides IRS spectra observed in our own programs 
(Program-ID 49, PI F. Low; Program-ID 82, PI G. Rieke; 
Program-ID 3624, PI R. Antonucci; Program-ID 20142, 
PI P. Ogle), we searched for archived spectra for objects 
in the three parent samples. Our sample is listed in Ta- 
ble [T] Fig. [1] compares the final three subsamples with 
their corres ponding parent samples. For the PG parent 
sample from [Schmidt fc GreenI (|1983l ). we exclude a non- 
quasar object PG 0119-1-229 and corr ect the redshift of 
PC 13 52-hOll to be 1.121 according to'B oroson fc GreenI 
(II992D . As shown in Fig.[l] we have included the whole 
PG parent sample at z <0.5. The quasar PG 2349-014 is 
not included in the original PG parent sample and this 
is why our PG subsample has one more object in the 
second redshift bin. For the 2MASS and 3CR subsam- 
ples at z <0.5 and z<1.0, respectively, about one third 
of the objects are included in this study. The subplots 
show that our 2MASS and 3CR subsamples are strongly 
biased toward high flux density at the wavelength where 
their parent samples are selected. 

2.2. Data Reduction 
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The spectra were obtained with the IRS using the stan- 
dard staring mode. The intermediate products of the 
Spitzer Science Center (SSC) pipehne S13.0.1, S13.2.0 
and S15. 3.0 were processed w ithin the SMART software 
package (|Higdon et al.ll2004D. For a detaile d description 
of the data reduction, see iShi et all (|20f)fih . Iffi^s et all 
dlOOe) and lBouwman et al. I (|2006f ). 

The sht widths of the short-low (SL) and long-low (LL) 
modules are 3'.'6 and lO'.'S, respectively. In order to mea- 
sure the star formation from the entire galaxy, we need 
to evaluate the extended IR emission outside of the IRS 
SL slit. The SL slit width is several hundreds of parcsecs 
for 3C 272.1 and 3C 274, 2-10 kpc for sixty-one objects 
(z <0.17) and >10kpc for the remaining objects. For 3C 
272.1, the MIPS image shows extended IR emission from 
the host galaxy and that this emission is thermal based 
on the extrapolation from radio data. The extended IR 
emiss ion of 3C 274 i s dominated by non-thermal emis- 
sion (jShi et al.ll2007l ) and is not related to star forma- 
tion. For objects with physical slit widths between ^2 
kpc and 10 kpc, a total of seventeen objects show ex- 
cess IR fluxes in the LL modules compared to the SL 
modules. However, the flux difference between the SL 
and LL modules can be caused by different slit-loss due 
to pointing errors, not necessarily by extended IR emis- 
sion outside the SL module slit. For 14 out of these 17 
objects, we obtained archived MIPS 24 /im images and 
measured the FWHMs of the radial brightness proflles. 
AH of them show FWHMs smaller than 3 pixels (the PSF 
has a FWHM of 2.4 pixels), implying that the excess IR 
fluxes in the LL modules are not due to extended IR 
emission from the host galaxies. For the remaining three 
objects without MIPS 24 ^m images, we use 2MASS K- 
band images and find that the excess flux of LL relative 
to SL for one object (PG 2304-1-042) may be due to ex- 
tended IR emission. For objects with slit widths larger 
than 10 kpc, we simply assume that the IRS slit con- 
tains all the IR emission from the galaxy and that the 
mismatch between the SL and LL spectra is due to vari- 
able slit-loss. Therefore, except for 3C 272.1, 3C274 and 
PG 2304-1-042, we rescale the SL spectra so that the SL 
and LL spectra have the same flux density at 14.2 /im. 

2.3. The Extraction of Aromatic Features 

The 7.7 /im feature resides at the blue end of the sil- 
icate feature. The level of contamination by the silicate 
feature on the aromatic flux measurement depends on 
several factors, including the strength of the silicate fea- 
ture, the shape of its blue wing and the shortest wave- 
length that the blue wing extends to. As shown in 
iHao et al.l (l2005l ) or our Figure [3l all these factors vary 
in different sources, resulting in deviations from the line 
profile for a normal galaxy interstellar medium. To ac- 
count for these variations, we fit the blue wing of the 
silicate feature with a Doppler profile: 



fx 



/Ao 



(A-Ao)2 + K): 



;exp{^{{\~\o)/aDf) (1) 



where Aq can be interpreted as the central wavelength 
of the silicate feature, and the combination of ao and 
aL control the shape of the blue wing and the starting 
wavelength where the silicate feature arises. The profile 
has no physical meaning and is adopted only for practical 
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Fig. 2. — Examples of the extraction of the 7.7 and 11.3 um 
aromatic features in the spectra with silicate emission, no silicate 
feature and silicate absorption, respectively. The dotted lines are 
the IRS spectra while the solid lines are the continua. The subplots 
show the Drudc profiles of the two features where the 11.3 um 
feature is fitted with two Drude profiles (dotted lines). 

purposes. As shown in Figure [U it can fit the 7.7 /im 
feature well. 

The procedure to extract the 7.7 /im aromatic feature 
is as follows. The spectra are first rebinned to a reso- 
lution of 0.1 /im to remove multiple points at the same 
wavelength, using the SMART software. The continua 
underlying the 7.7 /xm aromatic features and silicate fea- 
tures are defined as power laws over three spectral win- 
dows, 5.2-5.5 /J,m, 5.5-5.8 /im and 6.7-7.0 /xm. These 
spectral regions are selected to avoid the possible ice fea- 
ture at 6.0 /im and aromatic features at 6.2 /im. We 
then fit the continua-subtracted spectra simultaneously 
with two aromatic features at 7.7 and 8.6 /im and the 
silicate feature. The shapes of the aromatic features are 
assumed to be Drude profiles. Due to the low EW of 
aromatic features in AGNs, the FWHMs of the 7.7 and 
8.6 /im features are fixed at 0.6 and 0.3 /zm, respectively. 
The height of the 8.6 /im feature is also fixed to be one- 
third of that of the 7.7 /im feature. This relative height is 
similar to those of t wo average sp e ctra o f Hll-like nearby 
galaxies obtained bv lSmith et al.l (|2007f ). For the silicate 
feature, we fit only the blue wing with a Doppler profile. 
The starting wavelength of the spectral range for the fit 
is fixed at 6.5/im. We vary the red end from 9 to 12 
/im to have the best fit judged by visual inspection. For 
most of the sources, the measured aromatic fiux depends 
little on the selected red end wavelength. The feature is 




Rest-frame wovelength (um) 

Fig. 3. — IRS spectra of AGN with detected aromatic features. The soUd line is the derived continuum for the 7.7 /im and/or 11.3 /^m 
aromatic features. The subplots show the Drude profiles of the two features. 



considered detected if the height of the 7.7 /im feature is 
five times greater than the mean noise in the continuum. 

For the 11.3 fim feature, the sihcate feature behaves 
Hke a continuum and the slope of the underlying sili- 
cate profile varies smoothly across the aromatic feature. 
Therefore, we are able to determine the silicate profile 
simply with a quadratic interpolation. The 11.3 fiia fea- 



ture is fitted with two Drude profiles centered at 11.23 
and 11.33 with fixed FWHMs of 0.135 and 0.363 ^m, 
respectively. The combination of these two Drude pro- 
files fits well the 11.3 /zm features of nearb y galaxies, as 
demonstrated with high S/N IRS spectra bv lSmith et al.l 
(2007). After the spectrum is rebinned to a resolution of 
0.1 /im, the continuum (plus silicate feature) shape is 
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Fig. 3. — Continued. 



defined by using a quadratic interpolation over the four 
continuum spectral regions, 9.7-10.0, 10.0-10.3 /xm, 10.7- 
11.0 /um and 11.7-12.1 /xm. We then fix the continuum 

shape, the FWHM and the center wave-length of the two 
Drude profiles, but adjust the normalization of the con- 
tinuum and the strength of Drude profiles to fit the spec- 
tra in the range including the continuum and the feature 



(9.7-10.3 /im and 10.7-12.1 /im). The feature is consid- 
ered detected if the height of the combination of the two 
Drude profiles is five times greater than the mean noise in 

the eontinuTim. If the feature is not detected, the upper 
limit is calculated by assuming the same relative strength 
of the two Drude profiles as given by the fit and taking 
five times the mean noise for the total height of the two 
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Fig. 3. — Continued. 



Rest-frame wovelength (um) 



profiles. We visually inspected each detected feature and 
found that the 11.3 /im features of eleven objects may not 
be real due to larger noise around the feature relative to 
the mean noise in the continuum. For fifteen objects, the 
continuum was also fitted with an alternative quadratic 
interpolation, due to a large change in the slope of the 
silicate profile around the 11.3 /xm feature. However, the 



difference in the feature strength is smaller than a fac- 
tor of 1.5, showing that the continuum fitting procedure 
does not affect our results strongly. 

To test the robustness of our procedures against strong 
continua, power-law continua with differen t strengths 
are ad ded t o the star-forrn i ng te mplates from lDale et al.l 
(|200lh and lDale fc Heloul (l200l) . The 7.7 and 11.3 /^m 
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Fig. 3. — Continued. 



aromatic features are extracted using the above proce- 
dures and the flux variations are smaller than 1% for the 
EW range from the original value (~1/Ltm) to 0.01 um. 

2.4. Uncertainty of the Aromatic Flux 

We have evaluated each step in extracting the features 
to estimate the final uncertainty of the aromatic flux. 



To estimate the uncertainty due to the rebinned spectral 
resolution, the fluxes are re-measured with rebinned res- 
olutions from 0.08 to 0.12 /xm for features observed with 

SL module (resolution of ^0.1/im). For the objects at 
z>0.24, where the 11.3 /im feature is observed with LL 
module ( resolution of '^0.28/im), wc compare the mea- 
sured flux for rebinned resolutions ranging from 0.1 to 
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0.3 fiia. Comparing these measurements to the feature 
flux obtained at a rebinned resolution of 0.1 /im, we find 
that the differences are always below 10%. 

To estimate the uncertainty caused by the photon noise 
and the fit of the continuum and silicate feature, we pro- 
duce a noiseless spectrum for each detected aromatic fea- 
ture. The simulated noiseless spectrum for the 7.7 fim 
feature is the measured power-law continuum plus the 
measured Dopper profile of the silicate feature plus two 
Drude profiles of the measured 7.7 and 8.6 /xm features. 
The spectrum for 11.3 /xm is the quadratically interpo- 
lated continuum and silicate profile plus two measured 
Drude profiles. We then perturb this noiseless spectrum 
100 times to produce noisy spectra with mean S/N equal 
to the observed S/N. The aromatic features are then ex- 
tracted from these simulated spectra in the same way and 
the 1-cr uncertainty is obtained as the difference between 
the original flux and those from the simulated spectra. 
The uncertainty in this step is typically <15% for the 
11.3 /im feature and <30% for the 7.7 /xm feature. 

Due to the contamination by the silicate feature, we are 
unable to fit the 7.7 /im feature with multiple Drude pro- 
files. To compute the uncertainty in the assumed profile 
with a fixed FWHM for the 7.7 fim f eature, we ha v e used 
the code (PAHFIT.pro) written by ISmith et al.) (|2007l ) 
to measure accurate fluxes for the 7.7 /Ltm aromatic com- 
plexe s of the f o ur co mposite spectra of nearby galaxies 
in Sm ith et al.l ()2007f l. We then re-construct the 7.7 fim 
profile with the fitted parameters and measure the flux 



with a single Drude profile with a FWHM of 0.6 fim. The 
difference in fitted feature strengths is around 10%, which 
is adopted as the uncertainty due to the 7.7 fim aromatic 
profile. No uncertainty is applied for the assumed pro- 
file of the 11.3 fim feature. The above uncertainties are 
added quadratically to give the final error of the mea- 
sured aromatic flux. Table [1] lists the measured fluxes, 
uncertainties and EWs for both aromatic features. 

3. EXCITATION MECHANISM OF AROMATIC FEATURES 
IN AGNS 

As shown in § 1, the low EW of the aromatic features 
and the spatial extension of the aromatic emission in ac- 
tive galaxies suggest that these features are most likely 
predominantly excited by star formation. With the sig- 
nificant number of detections of aromatic features in this 
study, we can test this hypothesis. 

3.1. The Profile of Aromatic Features in AGN 

3.1.1. The Composite Spectra 

To study the profile of the aromatic features in AGN, 
we have produced the composite spectra for several 
groups of objects. The composite sp ectrum is computed 
follow ing the procedure described in IVanden Berk et all 
(|2001f ). All the observed spectra are shifted to the rest- 
frame and then rebinned to a common spectral resolu- 
tion (0.1/im) within SMART. After they are ordered by 
redshift, the first spectrum is rescaled randomly. The fol- 
lowing individual spectrum is rescaled to have the same 
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Fig. 4. — (a) The number of objects in each wavelength bin 
of the composite spectrum, (b) The arithmetic mean spectrum 
(the heavy solid line) of AGN with one of the 7.7 and 11.3 fim 
aromatic features detected and the fitted continuum (the light solid 
line), (c) The continuum-subtracted spectrum (the heavy solid 
line) superposed with the composite spectrum (the dotted line) of 
the Hll-like galaxies from lSmith et al.l l|2007l '). 

mean flux density in a common spectral region of the 
mean spectra of all lower rcdsliift spectra. The common 
spectral region is defined to be 5.0-6.0 fim where there 
is little influence from the silicate or aromatic features. 
The final composite spectru m is the arithmetic mean o f 
all rescaled spectra. Unlike in lVanden Berk et all (|2001h . 
we have not produced the median spectrum since the av- 
erage one shows much higher S/N. As implied by the 
compositing procedure, the aromatic features of individ- 
ual observed spectra with higher EW have larger weight 
in the feature of the final composite spectrum. 

The first arithmetic mean spectrum is the one of AGN 
with at least one of the 7.7 and 11.3 /im features detected. 
Fig. (Ha) plots the number of objects used in each wave- 
length bin. As shown in Fig. i the overall spectrum 
shows a power-law-like continuum with weak silicate fea- 
tures. We determined the continuum between 5.0 and 
10.0 /im using the procedure for extracting the 7.7 /xm 
feature but do not constrain the strength of the 8.6 /im 
aromatic feature. The continuum between 9.5 and 14.0 
/im is defined to be a quadratic interpolation over the 
mean flux densities of four spectral regions (10.0-10.3, 
10.8-11.0, 13.0-13.2, and 13.4-13.6 /im). As shown in 
Fig. life), broad features are present at 6.2 /im, 7.7 /im, 
8.6 /im, 11.3 fim and 12.8 /irn, similar to those in sta r 
forming galaxies fSee lLu et all l2003t [Smith et al.ll2007[ ). 



Fig. 5. — (a) The number of objects in each wavelength bin of the 
composite spectra of PG, 2MASS and 3CR AGN, respectively, (b) 
The arithmetic mean spectra and the fitted continua (the light solid 
lines), (c) The continuum-subtracted spectra of PG and 2MASS 
AGN, superposed with th e composite spectra (the dotted lines) of 
the Hll-like galaxies from lSmith et al.l l|2007l '). 

The dotted curve in Fig. ^c) shows the mean spec- 
trum of two compo site spectra of Hll-like galaxies from 
ISmith et al.l (|2007[ ) where the spectrum is shifted and 
rescaled to match the 11.3 /im feature of the AGN spec- 
trum. There is only a small discrepancy in the shapes 
and relative strengths of the aromatic features between 
AGN and Hll-like galaxies. A small amount of excess 
emission at 7.7 and 12.8 /im in the AGN spectrum is most 
Hkely due to [NeV]7.65/im and [NeII]12.8/im, respec- 
tively, as the excess emission shows a narrow FWHM. 
The result indicates the observed aromatic features in 
AGN resemble those in star-forming galaxies. The com- 
posite spectrum of AGN without either feature detected 
still does not show detectable aromatic features. 

Fig. El shows the arithmetic mean spectra for PG, 
2MASS and 3CR objects, respectively. The silicate emis- 
sion features are present in the PG spectrum while the 
2MASS and 3CR spectra have silicate absorptions. Aro- 
matic features are visible in the PG and 2MASS com- 
posite spectra, but not in the 3CR spectrum. As shown 
in Fig. ISI^c), the comparison to the Hll-like galaxies in- 
dicates the 11.3/7.7/im feature ratio (~0.30) of the PG 
spectrum is a little higher while the 2MASS spectrum 
presents a lower ratio (~0.22). However, they are within 
the one-(T range for star-forming galaxies as shown below. 

3.1.2. The Distribution of the Aromatic Feature Ratio 
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Fig. 6. — The ratio of the 11.3 /im aromatic flux to the 7.7 
flux. Th e uppe r pl ot is the ratio for n ormal spiral galaxies from 
ILu et al.l II2003I ') and lSmith et al.l II2007I) while the lower plot is for 
active galaxies in this paper. 



The above comparisons reveal that the shapes and 
relative strengths of the aromatic features of the AGN 
composite spectra are similar to those of Hll-like galax- 
ies. Fig. [S] compares the distribution of the 11.3/7.7/im 
aromatic rati os between AGN an d nor mal star- f ormin g 
galaxies from ISmith et al.) (I2007D and iLu et al.l (|2003f ). 



For the sample of ISmith et al. ~ |2007( ). we only include 



Hll-like galaxies but exclude a low-metallicity dwarf 
galaxy (HoII). No correction is applied to their aromatic 
fluxes, since they are obtained with multiple Drude pro- 
file fitting. The flux of the 7.7 an d 11.3 /xm aromatic fea- 
tures quoted in 'L u et al] (|2003( ) is the integrated value 
without continuum subtraction from 7.20 to 8.22 /im 
and from 10.86 to 11.40 /im, respectively. We correct 
their ratios by a factor of 1.08 to account for the differ- 
ence between their measured fluxes and the Drude-profile 
fluxes used in this paper. This factor is obtained based 
on the four compo site spectra of ne arby g alaxies from 
ISmith et al.1 (|2007D . In the iLu et all (|2003f) sample, one 
object is excluded since the integrated aromatic flux con- 
tains signiflcant hot dust emission. 

As shown in Fig. [6l the flux ratio of AGN with both 
features detected has a similar distribution to that of 
star-forming galaxies. The mean 11. 3/ 7. 7- aromatic ra- 
tio for the AGN is 0.27±0.1, compared with 0.2 8 ±0.11 
and 0.26ib0.07 for the spiral galaxies of iLu et all (|2003D 
and 'Smith et al.l (|2007t l. respectively. The Kolmogorov- 
Smirnov (K-S) test indicates a probability of 99% and 
40% that our AGN sample is the same as the s t ar-for ming 
galxies of iLu et al.l (|2003f ) and ISmith et all (|2007t l. re- 
spectively. 

Variations of the aromatic flux ratio have been ob- 
served among regions covering a wide range of physi- 
cal and chemi cal properties (e.g. Roclfse ma et al.lll996l : 
I Vermeil et al.|[2002.) . On the other hand, studies of the 
aromatic features in the same environment show that 



the flux r atio is insensitive to the intensity of t he radia- 
tion fleld (jUchida et al.ll2000l : IChan et al.ll2001h . Among 
different galaxies, there is no systematic difference in 
the aromatic flux rati o with th e inten sity of the radia- 
tion field, as seen by iLu et al] ()2003l ) where the spiral 
galaxies studied have total IR luminosity spanning from 
10^ to 10^^ Lq. This may arise because various aro- 
matic regions are averaged out over the entire galaxy. 
The similar distribution of the ratio between AGN and 
spiral galaxies as shown in Fig. (5) implies that the fea- 
tures observed in AGN are excited under conditions sim- 
ilar to those averaged over normal star forming galaxies. 
ISmith et all ()2007l ) have found that 20% of galaxies with 
low-luminosity active nuclei show a weak 7.7 fim feature 
relative to the strength of the 11.3 /im feature. The ori- 
gin of this deviation is not well understood. However, if 
it is the nuclear radiation that accounts for this peculiar 
ratio, the similar feature ratio between our sample and 
star-forming galaxies indicates the aromatic feature out- 
put in our sample is dominated by star formation, not 
by the active nuclei. For objects with only one detected 
feature, the distribution of the limits on Ft^t^jti/ Fii.^iim 
is still consistent with that of star-forming galaxies. 

3.2. The Global IR SED 

The global IR SED of AGN is affected by many fac- 
tors. However, if the aromatic feature originates from 
star-forming regions, the composite spectrum of the sub- 
sample with a higher fraction of aromatic emission in the 
mid-IR emission should show a higher fraction of far-IR 
emission. 

Fig. |7| compares the composite spectra from 5 to 200 
/zm for high-L(PAH)/L(MIR) and low-L(PAH)/L(MIR) 
objects, where i(MIR) is the total mid-IR luminosity 
between 5.0 and 6.0 /im and L(PAH) is the 11.3 /im 
aromatic luminosity or the 7.7 /ini aromatic luminosity 
multiplied by a factor of 0.27 for objects with only the 
7.7 //m feature detected. We define the dividing value 
of L(PAH)/L(MIR) for aU objects with MIPS 70 /im 
measurements so that the high and low-L(PAH) /L(MIR) 
subsamples have similar numbers of objects. The objects 
with upper limit measurements for the aromatic fiuxes 
are also included for the low-L(PAH)/L(MIR) subsample 
while only feature-detected objects are included for the 
high-L(PAH)/L(MIR) subsample. We have produced 
geometric mean composite s pectra, which conserve th e 
global continuum shape (See IVanden Berk et al.ir2001h . 
For each subsample, the IRS spectra are redshifted and 
rebinned to a common spectral resolution (0.1/im). The 
MIPS fluxes are K-corrected by assuming a=l and a=0.0 
{iu oc t/ ~°), respect i vely, based on the IR SEP of the 
AGN in ' Haas et all (l2003i) and iShi et all (l200l . Then 
each spectrum is normalized by the mean flux density in 
the wavelength range between 5.0 and 6.0 /im. The final 
composite spectrum is defined as (11" /a. i)^^" where A 
is the wavelength of a wavelength bin and n is the total 
number of spectra in this bin. 

Fig. iTja) plots the number of objects in each wave- 
length bin. As shown in Fig. |7l(b), given that the two 
composite spectra have similar weak silicate features, ob- 
scuration does not account for the difference in the shape 
of the SEDs. The high-L(PAH)/L(MIR) subsample has 
relatively larger IR emission toward wavelengths longer 
than 15 /xm. /(70/im)//(5-6/im) and /(160/im)//(5- 
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Fig. 7. — (a) The number of objects in each wavelength bin of the 
composite spectra of the high-L(PAH)/L(MIR) subsample (solid 
line plus filled circles) and the low-L(PAH)/L(MIR) subsample 
(dotted line plus open circles), where L(MIR) is the total mid-IR 
luminosity between 5.0 and 6.0 fira and L(PAIf) is the ll.Sfim aro- 
matic luminosity or the 7.7/xm aromatic luminosity multiplied by a 
factor of 0.27 for objects with only the 7.7 fim feature detected, (b) 
The geometric mean spectra of the two subsamples. (c) The spec- 
trum of high-L(PAH)/L(MIR) minus low-L(PAH)/L(MIR) objects 
superpos ed on the s t arbur st te mplate with L8_iooo um=2.0xlO-'^-^ 
Lq from lDale eTall ll200ll) and lDale fc Heloul | |200^ . 

6/im) are redder by a factor of 2.5 compared to the val- 
ues for the low-i(PAH)/L(MIR) subsample. The redder 
far-IR color is consistent with the star-formation origin 
of the aromatic features in these active galaxies. 

The spectrum of the high-L(PAH)/L(MIR) minus 
the low-L(PAH)/i(MIR) composite spectra is plotted 
in Fig. Ulc). We match th is residual spe c trum with 
star-forming templates from iDale fc Heloul (|2002l ) and 
find that the template with LiR(8-1000/xm)=2.0x 10"Lq 
presents the most consistent 70/160/xm color. After scal- 
ing this template to the 70 /xm photometry of the resid- 
ual spectrum, the subplot shows a good match for the 7.7 
and 11.3 fim aromatic features, although there is some 
discrepancy for the [NeII]12.8/im line. This match pro- 
vides further evidence for the star- formation origin of the 
aromatic features in these AGN. 

3.3. Molecular Gas 

Fig. [5] shows the mass of CO-derived molecular hy- 
drogen gas versus the aromatic-based star-forming IR 
(SFIR) luminosity (triangles). The aromatic-based SFIR 
luminosity is calculated in § 4. The mass of hydrogen gas 
is calculated using Mu:,^1.17AxlO^{Sco^V)Dl/ {1 + z), 
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Fig. 8. — The plot of the mass of CO-derived molecular hy- 
drogen gas versus the aromatic-based total IR luminosity (trian- 
gles) for AGN. Open and filled circles indicate weakly-interacting 
normal galaxies and strongly interacting normal galaxies from 
ISolomon &: Sag^ | |T988) . respectively. 

where 5*00^^^ is the CO flux in Jy km s^^ and £>l is 
the luminosity distance i n Mpc. The c ircles in Fig. [5] are 
the normal galaxies from ISolomon fc'S agc (1988), where 
open circles are for weakly-interacting normal galaxies 
and filled circles for strongly interacti ng ones. The to- 
tal IR luminosity LiR(8-1000/im) of the [Solomon fc Sagd 
()1988f ) sample is computed from IRAS four-band pho- 
tometry using the relation of Sa nders fc lylirab cl (199^ 
The difference between the relation of lSanders fc MirabeH 
([1996. ) and the star-forming templates used to derive the 
aromatic-based SFIR luminosity is typically less than 
5%. All physical parameters were corrected to our 
adopted cosmological model. Fig. [8] shows that the be- 
havior of the aromatic-based SFIR luminosity follows 
that of normal galaxies well. The relationship between 
the CO luminosity and SFIR luminosity is consistent 
with the star-formation excitation of the aromatic fea- 
ture in our AGN. 

As shown above, the profile of aromatic features, the 
global IR SED of AGN and the gas content in their host 
galaxies are all consistent with the predominantly star- 
formation excitation of the aromatic features in active 
galaxies. This conclusion confirms previous arguments 
based largely on s patially resolved s pectra of nearby ac- 
tive g aalxies (e.g . iCutri et "al]ll984l: [ Desert fc Dennefeic 



1988; [^ [l99a iLaurent et all l2000t iLe Floc'h'eTal 



2001 ). 



4. THE CONVERSION FACTOR FROM AROMATIC FLUX 
TO THE SFR 

Before proceeding with a quantitative study of the 
current star formation around AGN based on the 
measured flux of the aromatic features, we need to 
know how well the aromatic features trace the ongoing 
star- formation activity. For Galactic HII regions, the 
variation of PAH/far-IR(40-500/im) is up to two orders 
of magnitude from ul tra-compact to ext ended opti- 
cally visible examples ([Peeters et al.l [2003 ). However, 
integrated over the whole disk of spiral galaxi es, the 
aromatic features correlate well with Ha (Rouss eT et al.l 
I2001f) . This behavior may result from the galaxy-scale 
quantity averaging out the local physical properties in- 
volved in individual regions, such as the escape efficiency 
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of ion izing photons from HII regions (e.g. iRoussel et al.l 
1200 It ). The situation becomes comphcated in the 
circumnuclear regions where the EW of the observed 
aromatic feature is lo w, as in embedd e d HII regions 
(IRoussel et al.l 120011 : iHaas et al.l l2002t iPeeters et all 
I200J). The reason for this is unclear; it may be caused 
by oljscuration, PAH destruction, a decrease in ionizing 
photons as a result of the increasing compactness of 
the HII regions, or the additional mid-IR emission 
from highly embedded active nuclei. However, a direct 
attempt to correlate the aromatic feature to far-IR lumi- 
nosity for star-forming galaxies shows that the variation 
of PAH /far-IR is about a factor of 2-3 (IPeeters et al 



I2B0I iSpoon et al.ll200l iWu et al.ll200l . ISpoon' et al 
([2OOI obtained L(6.2/imPAH)/L(IR)=0.003±0.001 
from 70 normal and starburst galaxies. Taking a 

typical valu e of L(7.7/imPAH)/L(6.2/fmPAH)=3.5 

(jSmith et al.l l2007f ). this measurement is equivalent 
to L(7.7/imPAH)/L(IR)=0.01±0.0035. The aperture 
mismatch between the IR flux and the aromat i c flux 
contributes to a part of the scatter. iLutz et all (|2003l ) 
derived L(7.7/imPAH)/L(;8-1000Aim)=0.033±0.017 
(assuming a Drude profile with 0.6 /im FWHM for 
the 7.7 /J,m feature) from 10 starburst galaxies. This 
ratio allows for the aperture differences, although the 
two quantities are still not well matched. Ba s ed on 
IRS spectra of nearby galaxies, ISmith et al.l (|2007l ) 
employed a robust method of extracting aromatic 
features. The aperture-matched mean values with 
1-tT uncertainties of L(7.7/imPAH)/L(3-1100/j,m) and 
L(11.3/imPAH)/L(3-1100/im) are 0.052(1±40%) and 
0.012(1±30%), respectively, for 26 Hll-like normal 
galaxies excluding one dwarf galaxy (Ho II) with an 
extr eme low ratio prob ably caused by metallicity effects 
(See 'Smit h et al.ll2007t) . A part of the scatter in the ratio 
of L(PAH) / L(totlR) may arise from a general luminosity 
depen dence. As shown in Figure 3 of fSchweitzer et al.l 
(|2006f) . L(7.7/imPAH)/i/L,,(60Aim) decreases from 
0.06 for starburst galaxies at z^Liy(60/im)=1.5x 10^" 
Lq to 0.015 for starburst-dominated ULIRGs at 
vL^{60nm)=10^'^ Lq. 

To compute the luminosity-depende nt values, we have 
use d the star-form i ng tem plates from iDale et al.l ()2001l ) 
and IDale fc Heloul (|2002l ). Each SED template is opti- 
mized for a very narrow luminosity range (-^ ^ 0.1- 
0.4) where the luminosity is co nverted from the a in - 
dex using the relation given bv iMarcillac et al.l (|2006( ). 
Aromatic fluxes for all the templates are measured us- 
ing the same procedures as for AGN. As demonstrated 
in § 2.3, the aromatic fluxes obtained by our proce- 
dure do not change with the EW, implying that there 
is no systematic difference in the measurements of the 
aromatic fluxes between the star-forming templates and 
AGN. The conversion factor for the 7.7 /im feature varies 
from 0.041 at a SFIR luminosity of 10^ Lq to 0.0095 
at a luminosity of 3.3x10^^^ Lq and the 11.3 /im fea- 
ture varies from 0.012 to 0.004 over the same luminosity 
range. These values agree well with the observational 
ones. To derive the conversion factor for each object, 
we adopt the template that gives the closest aromatic 
flux at the redshift of this object. The uncertainties 
are assumed to be the observed ones (40% and 30% for 
L(7.7/imPAH)/L(8-1000/im) and L(11.3/imPAH)/L(8- 



1000/im), respectively), although there is only a 10% dif- 
ference between conversion factors for SED templates in 
two adjacent luminosity ranges. The final uncertainty 
of the aromatic-derived SFIR luminosity includes that of 
the conversion factor and the measurement uncertainty 
of the aromatic flux. If this final uncertainty is larger 
than the measured aromatic flux, the 3(T upper limit is 
adopted. Table [1] lists the SFIR luminosity calculated in 
the above way. For objects with both features detected, 
we adopted the value from the 11.3 /im feature since it 
generally has smaller uncertainty. The value from the 
detected feature is listed if only one feature is detected. 
For objects with neither feature detected, the lower value 
for the two upper limits is listed. 

As discussed in § 2.2, PG 2304-f 042 and 3C 272.1 have 
thermal IR emission outside the IRS slit. This extended 
emission is converted to the total IR luminosity by multi- 
plying by a factor of 12.0 based on the sta r-forming tem- 
te with L7ij(8-1000/im)=10" Lq fro m [Dalefc Heloul 
03), and is close to the observed value (jCharv fc Elbad 
'I). 

Non-star-formation sources, such as planetary nebulae 
and diffuse stellar radiation, can excite low-level IR emis- 
sion and aromatic features. Aromatic feat ures have been 
observed in a fraction of elliptical galaxies (jBressan et al.l 
2006) and some of them may originate from star forma- 
tion regions while others may be excited by an old stellar 
population. In five n ormal elliptical galaxies observed by 
iKaneda et al.l (|2005f l. the 11.3 /im aromatic luminosity 
is between 10^ and 8x10^ Lq (the possible problem in 
this work with stellar ligh t subtraction should not af- 
fect the 11.3 /im flux muchi iWregman et al.ll2006D . To be 
sure we are measuring recent star formation, we adopt a 
limiting aromatic luminosity of 3x10'' Lq above which 
the old stellar population contribution should be smaller 
than 25%. The corresponding aromatic-derived total IR 
luminosity at this limit is 3x10^ Lq. Therefore, a total 
of twenty-two objects including eight PAH-detected ones 
are excluded. 

5. ORIGIN OF THE FAR-IR EMISSION OF AGN 

Fig. ini shows the star-formation contribution to the 
MIPS rest-frame 24, 70 and 160 /im emission versus the 
integrated mid-IR luminosity between 5.0 and 6.0 /tm. 
The IRAS or ISO 25 /im fluxes are plotted for objects 
without MIPS 24 /im flux measurements. For objects 
without MIPS 70 fim flux measurements, we estimate one 
by interpolating between the detected IRAS or ISO 60 
and 100 /tm fluxes. The MIPS fluxes are K-corrected by 
assuming a=l for 24 and 70 fim photometry, and a=0.0 
for 160 /im photometry (f^ cx i^~°), based on the IR SED 
of AGN in IHaas et all (120031 ) and lShi et all (I2005tl . The 
total PAH-derived SFIR luminosities are converted to 
the star-formation emission at the three MIPS bands us- 
ing the luminosity-dependent conversi on factors derived 
from the star-form i ng te mplates from lDale et al.l (|2001f ) 
and Dale fc Heloul (|2003 ). 

At 24 /im. Fig. |5| indicates most of the objects are 
dominated by AGN emission. At 70 and 160 /tm, the 
far-IR emission of an individual AGN can be dominated 
by either AGN power or star formation. To quantify the 
star-formation fraction at the three MIPS bands and its 
possible dependence on th e AGN lumin osity, we have em- 
ployed the code written bv lKellvl (|2007l ) that incorporates 
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Fig. 9. — The star-formation fraction at 24, 70 and 160 
versus the mid-IR (5-6/im) luminosity for the PG, 3C and 2MASS 
objects, respectively (see the color version online). 

the upperlimit measurements. As listed in Table [2l the 
average star-formation fractions for the whole sample at 
MIPS 24, 70 and 160 fim are 4%, 26% and 28%, respec- 
tively, at the median mid-lR luminosity (2.6x10^*^ Lq) of 
the sample. As indicated by Tabled these ratios depend 
on luminosity, with a lower relative star-formation con- 
tribution at higher AGN mid-lR luminosity. The diverse 
nature of far-lR emission is consistent with the large scat- 
ter of the correlation betw een the far-IR emission and 
AGN power indicators [e.g.lShi et al.ll2005l : IClearv et all 
120061 iTadhunter et al.l 120071). There will also be some 
scatter due to the range of redshifts. However, since the 
redshifts of our PG and 2MASS samples are similar and 
modest, the effect should be small. 

TableHalso includes the resuh for PG and 2MASS ob- 
jects at the MIPS 24 and 70 /im bands, where there are 
enough detected data points. The average star-formation 
contributions at MIPS 70 ^m for PG and 2MASS are 
24% and 51% at median mid-lR luminosities of 3.0x10^° 
Lq and 3.5x10^° Lq, respectively. The fraction for 
the PG quasars is lowe r than that (>30%) obtained by 
iSchweitzer et al.l ()2006l ). who also employ the aromatic 
feature to evaluate the role of star formation. Contri- 
butions to the discrepancy include a difference in the 
conversion factors from the aromatic fluxes to the SFIR 
fluxes and the relatively large uncertainties in their 7.7 
fim fluxes caused by silicate features, whereas our result 
is mainly based on 11.3 /im features. 

Compared to the whole sample, PG objects show 



Fig. 10. — The star-formation fraction at 24, 70 and 160 fim 
versus the ratio of mid-IR (5-6 ^tm) luminosity and the Eddington 
luminosity for PG quasars. The solid line is the regression line and 
the two dotted lines are 2cr confidence bounds. 

relatively stronger luminosity-dependence of the star- 
formation fractions at 24 and 70 /im, with decreasing 
fractions at higher mid-IR luminosities. However, the 
2MASS objects do not have such a relation and most of 
the 3CR results are upperlimits. Thus the relation for 
the whole sample is mainly produced by the PG sample. 
As shown in Fig.fTUland Table[31 the star-formation frac- 
tions for the PG objects also decrease as the ratios of the 
mid-IR continuum luminosities and the Eddington lumi- 
nosities decrease, where the blackhole masses of PG ob- 
jects a£e_obtamed^wn IVest ergaard & Peterson (200(|) 
and iKaspi et al.l (j2000D . The anti-correlations indicate 
these two relations are not caused by the selection ef- 
fect that the detectable aromatic features in objects with 
higher mid-IR continuum emissions have larger fluxes. 

6. STAR-FORMING IR LUMINOSITY FUNCTION OF 
QUASAR HOST GALAXIES 

6.1. Methodology 

The main challenge in deducing the SFIR luminosity 
function (LF) for our sample is that the flux limit of the 
aromatic feature is not well defined and many objects 
have only upper limits in these measurements. Therefore, 
we obtained the SFIR LF by converting the well-defined 
LF at other wa velengths using the fractional bivariate 
LF (jElvis et al.i .l978). The formula can be written as 



MsFIR 



(2) 
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where $j 



is the SFIR LF and ^m, is the LF at A- 



A^sFiR 1° '^'^^ or in i^r ana 'i'Mx 
band where each parent sample is selected (radio for 3CR 
objects, B-band for PG objects and if-band for 2MASS 
objects). The fractional bivariate LF F{M\, MsFm) in- 
dicates the fraction of objects with magnitude M\ at 
A-band having SFIR luminosity of Msfir- We calcu- 
late f (Ma,Msfir)AMaAMsfir by dividing the num- 
ber ni of objects with A-band magnitude in the inter- 
val Mx±AMx/2 and the SFIR luminosity in the interval 
Msfir±AMsfir/2 by the number n2 of objects with A- 
band magnitude in the interval M\±AAI\/2 that could 
have had detected aromatic features if they had SFIR lu- 
minosities of Msfir. ni is the observed number. Any ob- 
ject with A-band magnitude in the interval M\±AM\/2 
will be counted into n2, if it has a limiting SFIR lumi- 
nosity lower than Msfir. The limiting SFIR luminosity 
is defined as the minimum star formation rate to detect 
the aromatic feature (see § 2.3) plus any extended IR 
emission. 

For PG quasars , $/v/f, is the i3-band LF at 0.0< z <0.5 
from Table 9 of Schmidt fc Greeiil ()1983D , where the me- 
dian redshift of 0.25 is adopted to convert the appar- 
ent magnitude to the absolute magn itude and the K- 
correc tion is the same as described in ISchmidt fc GreenI 
(|1983l ). This B-hand LF has data coverage for Mb from 
-21.4 mag to -26.4 mag. A double-expon ential model 
(for the formula, see iLe Floc'h et~al] l2005| ) fits the B- 
band LF well and it is used to derive the for any 

given Mb between -21.0 and -26.5 for our PG subsam- 
ple. The SFIR luminosity of this PG subsample spans 
the range from 3.1x10^ to 2.4x10^^ Lq. To construct 
the fractional bivariate LF {F{AIb, Msfir)), the entire 
ranges of Mb and SFIR luminosity are each divided into 
four intervals. The final fractional bivariate LF {F{Mb, 
Msfir)) along with Poissonian uncertainties is listed in 
Table II 

For 2MASS o bjects, the LF at K band from 
iCutri et al] (|2001h is adopted as ■ A two-exponential 
model does not fit the data well and thus we interpolate 
the measured data points to get the space density at a 
given /iT-band magnitude. Table \5\ lists the final frac- 
tional bivariate LF {F{Mb, Msfir)) for 2MASS objects. 

For 3CR objects $Ma is the LF at 151 MHz from 
IWillott et al.l ()2001h . where the LF is obtained based on 
the 3CRR, 6CE and 7CRS samples. We use the ana- 
lytic LF of model C for a cosmological model of 17^=0, 
JIa^O and iJo=50 kms~^Mpc^^, because the LF for 
this cosmological model is close to that f or our cosmo- 
logica l model except for the Hq value. (jWillott et all 
l2001t l. We convert to our cosmological model by setting 
$i(Xi,z)dI^ = ^2(L2,z)dV2 jPeacock 198S). The radio 
luminosity at 151 A4Hz for our 3CR subsample is calcu- 
lated and K-corrected u sing the flux density and spectral 
index at 178 MHz from lSpinrad et all (| 1985( 1. Again, we 
limit our 3CR subsample to the redshift range between 
0.0 and 0.5 to match the PG and 2MASS redshift ranges. 
The final fractional bivariate LF {F{Mi]^imhz, Msfir)) 
with Poissonian uncertainties is listed in Table [6l 



6.2. Star-forming IR Luminosity Function of Active 
Galaxies 

6.2.1. Comparison to Field Galaxies 
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Fig. 11. — Star-forming infrared luminosity functions for the PG, 
2MASS and 3CR AGN. The dotted Une is the re-normaUze d lumi- 
nosity function of local field galaxies from ILe Floc'h et al.l (2005). 

The most important result from the fractional bivari- 
ate LFs in Table H Table E] and Table M is that objects 
with a large range of nuclear activity have a non-zero 
probability of having a high SFIR luminosity. The form 
of the fractional bivariate LF implies that SFIR LF of 
AGN host galaxies is much flatter than the LF of the 
AGN themselves. 

Fig. [n] shows the results for the SFIR LF for the PG, 
2MASS and 3CR subsamples. Each subsample has a 
brightness limit at the wavelength where it is selected. 
We set Mb <-21 for the PG subsample and Mk <- 
25.5 for the 2MASS subsample and Lisimhz >2x10'^* 
W Hz-i Sr-i for the 3CR subsample. The dotted line 
shows the re-normalize d IR LF of local field galaxies from 
ILe Floc'h et"an (pOOSi ) based on the IRAS and ISO re- 
sults; it agrees we h with p revious studies of the IR LF of 
field galaxies (See lRiekefc Lebofskvlll986HSanders et al.l 
[2003l) . In Fig. [TTl the SFIR LFs of the three subsam- 
ples are much flatter than the re-normalized LFs of field 
galaxies. 

We need to be sure that the flatter LFs are not just 
a result of the difficulty in measuring the SFR around 
a bright quasar. We first use Monte-Carlo simulations 
to test the robustness of the methodology used to derive 
the SFIR LF of AGNs. The following steps are taken 
to construct a sample that mimics the PG subsample: 
(1) a total of A^ofcj (>10000) objects is created over the 
redshift range between O.OOl and 0.5; (2) the comoving 
number density is constant over the redshift range; (3) a 
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i3-band luminosity within the range of the PG subsam- 
ple is assigned to each object randomly but the relative 
distribution is the same as the PG B-band LF; (4) simi- 
larly, each object has a randomly assigned IR luminosity 
with relative distribution defined by the SFIR LF of the 
PG subsample. In this case, the IR flux is not correlated 
with the i3-band flux; (5) a well-defined fiux limit is ap- 
plied in the i?-band while the IR fiux limit is randomly 
distributed over the whole range of the SFIR fiuxes of 
the PG subsample. After producing the above set of 
objects, the fractional bivariate LF is calculated based 
on those objects detected in the B-band. The final de- 
rived SFIR LF using the fractional bivariate LF follows 
the pre-defined SFIR LF within the Poission noise. The 
same result is obtained for the simulation in which the 
IR fiux is tightly correlated with the S-band fiux. 

Unlike the PG subsample, which is complete, the 
2MASS and SCR subsamples only contain one-third of 
their parent samples at z <0.5. To test for the effects 
of the sample incompleteness, we use only one third of 
the objects brighter than the _B-band limiting fiux cre- 
ated in the above simulations, with these objects having 
the brightest apparent i?-band magnitude. Again, the 
derived SFIR LF is consistent with pre-defined SFIR LF 
within the Poission noise. We also test using the one- 
third of the objects with the most luminous absolute B- 
band luminosity. The shape of the derived IR LF does 
not change but the normalization becomes smaller. The 
same result is obtained if the i?-band flux correlates with 
the IR fiux. Thus, for all three subsamples, the Monte 
Carlo code demonstrates the robustness of our method- 
ology to derive the SFIR LF of AGNs 

Because AGNs have strong mid-IR continua, aromatic 
features are detected only in host galaxies with intense 
star formation. We can now use the Monte-Carlo sim- 
ulation to demonstrate that this selection effect cannot 
account for the large difference in the SFIR LF between 
the field galaxy and PG quasars. In the simulation, we 
assume the SFIR LF of PG quasars actually follows that 
of field galaxies. For each PG object, we obtain the IR LF 
of field galaxies at the redshift of this o bject by assum- 
ing th at the local field galaxy IR LF from lLe Floc'h et al.l 
pOOSi) evolves with redshift as L*{z) = L*(0)(1 + zf-'^ 
and <I>*(z) = <i>*(0)(H-2:)° ''. We then randomly assign a 
SFIR luminosity to this PG object with a relative proba- 
bility that follows the LF of field galaxies at this redshift. 
The range of the simulated SFIR luminosities is from 
3.1x10^ to 2.4x10^^ Lq, consistent with the observed 
range for the PG quasars. Also, we assume that the to- 
tal probability in this luminosity range is equal to 1. In 
this case, all simulated IR luminosities are above the low 
luminosity cut (3x10^ L©), and thus bias the results to- 
ward the high luminosity end. Combining the simulated 
SFIR luminosity and the observed uncertainty or upper 
limit for each PG object, we can calculate the detection 
fraction for the aromatic features. After one thousand 
simulations, we find (despite the bias toward high lumi- 
nosity) that the detection fraction is only (28±3)%, much 
smaller than the observed value (48%). This large dif- 
ference indicates that our result is not simply due to the 
selection toward high levels of SFR caused by the AGN 
emission. 

We further measure the probability of producing the 
observed curvature of the SFIR LF if the PG quasar sam- 
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Fig. 12. — Star-forming infrared luminosity functions of PG 
quasars as a function of quasar brightness. The dashed line is 
the re-normalized luminosity f unction o f star formation in CfA 
Seyfert 1 galaxies from Maiolin o et al.l II1995I ). The dotted line 
is the re-normal i zed lu minosity function of local field galaxies from 
ILe Floc'h eTal] l|200i). The solid lines are Schechtcr-function fits 
to the two PG subsamples. 

pie actually has a field galaxy SFIR LF. In each simu- 
lation, all PG objects are assigned randomly SFIR lu- 
minosities as described above. Using the simulated lu- 
minosities and the observed uncertainties or upperlimits, 
a SFIR LF is constructed using the same procedure in- 
cluding the number of luminosity bins as the observed 
LF. All data points produced in a total of ten thousand 
simulations are rebinned to the same bins as for the ob- 
served PG LF. In four luminosity bins, the fractions of 
simulated non-zero number densities are 100%, 100%, 
64% and 6% from low to high luminosity. All simulated 
number densities are then rescaled by a factor to match 
the composite number density in the first luminosity bin 
to the observed one. This composite number density is 
assumed to be the median value of all simulated number 
densities (including zero value) in the first bin, indicating 
a probability of 50%. We then calculate the probability 
for an observed luminosity bin as the fraction of sim- 
ulated number densities larger than the lower 1-sigma 
bound of the observed number density in this bin. The 
probability in each bin from low to high luminosity is 
99.0%, 1.0%, 2.5% and 4.0%, respectively. This result 
provides further evidence that the flatter SFLF of the 
PG quasars is robust against selection effects. 

6.2.2. Dependence on AGN Luminosity 

Fig. [12] shows the SFIR LF of PG quasars as a func- 
tion of the i3-band luminosity. The two solid lines are 
Schechter-function fits for PG quasars at Mb <-21 and 
Mb <-23, respectively. The fitting parameters arc given 
in Table El There is a trend that the SFIR LF of PG 
quasars becomes flatter for the brighter PG objects. We 
suggest that the higher SFR for brighter PG quasars 
is not a selection effect because the S-band luminosity 
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of normal infrared galaxies is not well correlated with 
IR luminosity and LIR Gs rarely have Mb < -23 (See 
iRieke fc Lebofskvll 19861) . The trend seen in Fig. [His not 
likely to be due to evolution with redshift, as the mean 
redshifts for the faint and bright subsamples are nearly 
the same from faint to bright, 0.18±0.30 and 0.24±0.11 
respectively. 

In Fig. [121 the dashed line is the LF o f extended star 
forma tion in CfA Seyfert 1 galaxies from iMaiolino et al.l 
(|1995( ). The extended IR emission of Seyfert galaxies 
was obtained by subtracting the nuclear emission from 
IRAS 12 /zm photometry (See IMaiolino et al.l[l99l . We 
converted the 10 /im luminosity to th e total IR luminos - 
ity using the IR SEP te mplate from [Dale et all (|2001[ ) 
and [Dale fc Heloul (|200 21 . Similarly to converting the 
aromatic flirx to the total IR luminosity, the conversion 
factor from 10 /zm to the total IR luminosity depends 
on the total IR luminosity. The omis sion of nuclear 
star fo rmation (within 2") in the study of lMaiolino et al.l 
(|1993) may affect the LF of total star formation in their 
Seyfert galaxies. However, if nuclear star formation is 
correlated with t he ext ended star formation as found by 
[Buchanan et al .l (|2006') , the shape of the LF for the total 
star formation in Seyfert galaxies should not change. As 
shown in Fig. [121 the SFIR LF of Seyfert 1 galaxies is 
steeper than the LF of PG quasars. There is also a sug- 
gestion that the LF for the lower-luminosity PG quasars 
is steeper than for the higher-luminosity ones. Seyfert 
galaxies have a higher SFR and flatter LF on averag e 
than field galaxies (see Fig.[T^and'Maiol ino et aLlilQgSO . 
It appears that star formation is correlated with the level 
of nuclear activity over the full range from normal galax- 
ies to quasars. 

To test the trend of the SFIR LF of active galax- 
ies as a function of AGN luminosity, we extended the 
Monte-Carlo simulations described in § 6.2.1 to test the 
difference between PG quasars with Mb <-21 and PG 
quasars with Mb <-23. In this simulation, we assume 
that the SFIR LF of PG quasars with Mb <-23 actu- 
ally follows that of PG quasars with Mb <-21. For a 
PG quasar with Mb <-23, we obtain the SFIR LF of 
PG quasars with Mb <-21 at the redshift of this object 
by assuming the SFIR LF of PG quasars at Mb <-21 
evolving with redshift as L*{z) — L*{zi){j^Y-'^ and 

(^*[z) = $*(zi)(p^)° '', where zi is the mean redshift 
(0.2) of PG quasars with Mb <-21. Based on this LF, 
a random SFIR luminosity is assigned to a PG quasar 
with Mb <-23. The luminosity range is between 3.1 x 10^ 
and 2.4x10^^, consistent with the observed range for PG 
quasars with Mb <-21. The total probability in this lu- 
minosity range is equal to 1. Using the observed uncer- 
tainties or upper limits, we predict the detection fraction 
of the aromatic feature for PG quasars at Mb <-23 of 
17±5%, smaller than the observed fraction of 28%. This 
result supports our conclusion that the SFIR Imninosity 
increases with increasing AGN luminosity. 

6.2.3. Comparison Between Different Subsamples 

As shown in Fig. [Til the behavior of star formation is 
different around AGN selected by different techniques. 
Since the SFIR LF of AGN host galaxies is a function of 
AGN luminosity as found in the last section, the effect of 
the nuclear brightness needs to be removed. The 2MASS 
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Fig. 13. — Cumulative fraction luminosity functions F(>L) = 
^1^=1,0 ^^^^ ^'^^ objects versus 2MASS objects (upper plot) 
and the PG objects versus 3CR objects (lower plot), where f(L) is 
the fractional luminosity function (See text). 



iiT-band photometry for all PG objects was obtained from 
the 2MASS Point Source Catalog. We calculated B-K 
for all PG objects and found that < B-K >=3.0±0.6 
and is not a function of absolute if -band magnitude. 
All PG objects with Mb <-22.5 are selected to form a 
comparison sample for the 2MASS objects with Mk <- 
25.5. For the SCR subsample, it is difficult to select a 
PG sample with the same level of nuclear activity. This 
is because PG objects are selected by thermal emission 
while 3CR objects are selected because of their non- 
thermal emission and there is no good correlation be- 
tween the radio emission and the thermal mid-IR emis- 
sion (Ogle et al. 2006). Instead, we compare the whole 
PG subsample at Mb <-21 to the whole 3CR subsample 
at Li5iMHz >2xl024 W Hz-i Sr'^. 

Fig. [13] shows the cumulative fractional luminosity 
function F(>L) = T,T=Lo for PG versus 2MASS 
and PG versus 3CR. To avoid biases due to evolu- 
tion, the comparison includes objects with z <0.5. The 
fractional luminosity function f(L) is defined similarly 
to the fractional biv ariate LF (See I Elvis et al] 119781 : 
iGolombek et al.lll988l ). As shown in Fig. [131 there is an 
apparent sequence in terms of the level of SFR that pro- 
gresses from 3CR to PG to 2MASS objects that generally 
show the highest SFRs. The median star-forming IR lu- 
minosities of 3CR, PG and 2MASS objects are 6x10^, 
3.0x10^° and IxlO" Lg, respectively. Different AGN 
selection techniques appear to identify objects with dif- 
ferent levels of star forming activity in their host galaxies. 
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6.3. Implications for Nuclear Activity 

The flatter SFIR LF of AGN host galaxies indicates en- 
hanced star-forming activity relative to local field galax- 
ies. Previous studies illustrate the presence of significant 
post-starburst stellar populations in quasar host galaxies. 
For example, the optical and near-IR broadband SEDs of 
AGN indicate the presence of young stellar populations 
with an age of about a Gyr in the host gal axies, inde- 
pendent of morphological tvpe d Jahnke et al.ll2004) , con- 



sistent with previous studies (iKotilainen fc Wardl UJO^ 
Schade et al.ll20qa iRonnback et al.lll996D . In addition. 



K auffman n et al. 1 12001 " bund a trend of younger mean 
stellar population for higher-luminosity AGN based on a 
very large sample. None of these studies found evidence 
for intense on-g oing massive star fo rmation, except for a 
few objects fsee lJahnke et al.ll2004D . We emphasize that 
the techniques employed in the above studies are unable 
to detect OB stars or suffer from strong degeneracy be- 
tween the current star-formation and the star-formation 
history. Therefore, these studies do not contradict our 
result. Searches for massive star formation through UV 
spectroscopy or spatially-resolved observations for star- 
formation tracers (such as recombination lines and IR 
emission) indicate the prese nce of massive star formation 
in Se yfert galaxies fMaiolin o et al.lll995HHeckman et al.l 
Il997f) and in quasars (,Cresci et al.ll2004f) . All of these 
studies focus on the central region of the galaxy, imply- 
ing that the star formation in quasars is circumnuclear. 
This is consistent with the lack of spectroscopic evidence 
for on-going star formation at distances from the nuclei 
of ~15 kpc (|Nolan et al.ll200lD . 

The flatter SFIR LF of AGN host galaxies relative to 
field galaxies also implies that nuclear activity tends to 
be triggered in galaxies with enhanced star formation. 
Based on Fig. [Til we can calculate the probability of trig- 
gering a PG quasar in field galaxies at a given SFR; for 
example, the probability of triggering nuclear activity at 
isF/i?,= 1-25x10^^ Lq is a factor of 50 higher than that 
at L5F//s=lxlO^° Lq. This indicates an environment 
with intense star formation offers preferential conditions 
for nuclear activity, such as th e abundant inflowing ma- 
terial driven by star formation (jGranato et al.|[2004 ). On 
the other hand, it implies that over much of the life of an 
AGN, its feedback does not quench the star formation, 
but instead may enhance the host galaxy star fo rma- 
tion a s demonstrated in some numerical simulations (jSilkl 
|2005( ). Our result that more luminous AGNs are more 
likely to reside in host galaxies with more intense star 
formation provides further evidence that feedback from 
the two physical processes (star formation and nuclear 
activity) can enhance both processes. Numerical simula- 
tions have predicted the evolution of the S FR and SMBH 
accretion rate a long the merging process (jGranato et al.l 
12004 ISpringeT ct al. 2005). They conclude that the evo- 
lution of star formation almost follows the SMBH accre- 
tion rate, although the former starts to decline a little 
earlier. A more quantative and careful comparison be- 
tween the simulations and our observations will improve 
our understanding of when and how feedback plays a role 
in galaxy evolution and SMBH growth. 

Although PG, 2MASS and SCR AGN have flatter 
SFIR LFs compared to field galaxies, they show dif- 
ferences in the distribution of SFRs, as indicated by 



the cumulative fractional LFs in Fig. [T3l Fig. [8] shows 
that the SFR of AGN host galaxies correlates with the 
amount of molecular gas in the host galaxy, which sug- 
gests that different AGN selection methods prefer host 
galaxies with different levels of gas reservoir. It is 
interesting that PG and 2MASS quasars have differ- 
ent levels of SFR. Both samples are selected through 
thermal emission. There is no obscuration along the 
line of sight for PG objects while the red IR-optical 
color of 2MASS objects is attributed to the obscu- 
ration of nuclear radiation by dust in the circumnu- 
clear regions or h ost galaxies (e.g. ISmith et al.l 120021 : 
iMarble et ani2003l|. Ac cording to the AGN unification 
model ( Antonuccil Il993[ ) . 2MASS objects are reddened 
counterparts of PG objects. The different levels of star 
formation in 2MASS and PG objects suggest that star 
formation affects our view of the AGN phenomenon, 
which is not expected under the unification model. This 
is not a selection effect that 2MASS objects need to 
have a larger SFR to have comparable the i^-band lu- 
minosity to PG quasars, as if -band fluxes in 2MASS 
objects are dominated by hot dust or starlight, not by 
star formation. A similar correlation has been observed 
in Seyfert galaxies, that Seyfert 2 objec ts have larger 
star formation rates than Seyfert Is (e.g. 'Edelso n et al.l 
Il987t [Maiolino et al. 1995). Observations and numeri- 
cal simulations show that the feedback produced by nu- 
clear star formation can heat the circu mnucl ear mate- 
rial and thus increase its scale height (iMaiol ino et aD 
'Ohsuga & Umemura"199i; IWada fc NOT^anii200l 
Watabc & Umcmura 2005). Such behavior could pro- 
duce the link between star formation activity and AGN 
properties. 

7. CONCLUSIONS 

We present Spitzer IRS observations of three AGN 
samples including PG quasars, 2MASS quasars and SCR 
radio-loud AGNs. The PG sample includes aU PG 
quasars at z<0.5 while one third of the 2MASS and SCR 
parent samples are used in this study. The main results 
are the following: 

1. The aromatic features at 7.7 and 11.3 /Ltm are de- 
tected against the strong mid-IR continuum of the AGN. 
The excitation mechanism for the aromatic features is 
predominantly star formation. 

2. The contribution of star formation to the far-IR 
emission of individual AGN is diverse; the average contri- 
bution is around 25% at 70 and 160 ^m. For the PG ob- 
jects, this contribution shows anti-correlations with the 
mid-IR luminosity and the ratio of the mid-IR continuum 
and the Eddington luminosity. 

S. The star- forming IR luminosity functions of AGNs 
are flatter than that of held galaxies, implying the feed- 
back from star formation and nuclear activity can en- 
hance both processes. 

4. The star-forming IR luminosity function of AGNs 
is correlated with the level of nuclear activity over the 
whole range from normal galaxies to bright quasars, with 
higher star formation rates for more intense nuclear ac- 
tivity. The 2MASS, PG and SCR AGNs have distribu- 
tions of star formation that follow the progression (from 
high to low SFR) of 2MASS-PG-SCR, implying that var- 
ious AGN survey techniques select host galaxies with dif- 
ferent levels of star forming activity. 
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Shi et al. 



TABLE 1 

AGN WITH ASSOCIATED PHYSICAL PARAMETERS 



source 


Redshift 


F(7.7/im) 


EW(7.7;im) 


F(11.3/im) 


EW(11.3;im) 


^SFIR 


^5— Gum 


ScoAF 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


PG0003+158 


0.450 






< 0.13 






1.2x 10 






PG0003+199 


0.025 


< 1.39 




0.29±0.06 


0.01 


(8.8±4.03)xl0 


3.3x10"'^ 






PG0007+106 


0.089 


< 1.39 




0.51±0.06 


0.03 


(3.2±1.32)xl0-'-'' 


1.9x10-"' 


< 3.00 


1 


PG0026+129 


0.142 


< 0.36 




< 0.12 




<4.5xl0^'' 


3.0x10-"' 






PG0043+039 


0.385 


< 0.36 




< 0.08 




<5.1xl0^^ 


l.Ox 10-'-' 






PG0049+171 


0.064 


< 0.50 




< 0.05 




<3.1xl0 


2.1x10"-' 






PG0050+124 


0.061 


8.28±5.61 


0.05 


2.77±0.25 


0.02 


(9.3±3.82)xl0-'-'-' 


4.3x10-"-' 


18.00 


2 


PG0052+251 


0.155 


< 1.74 




0.55±0.12 


0.05 


(1.3±0.62)xl0-'^-'^ 


3.2x10^" 


2.00 


3 


PG0157+001 


0.163 


6.71±2.44 


0.25 


2.44±0.16 


0.09 


(8.9±3.61)xl0^^ 


5.7x10-"' 


8.10 




PG0804+761 


0.100 


< 1.75 




< 0.19 




<3.8xl0^'' 


5.6x10^" 


2.00 


2 


PG0838+770 


0.131 


1.46±0.60 


0.17 


< 0.23 




(6.1±3.09)xl0^'' 


1.2x10^'^ 


3.40 


1 


PG0844+349 


0.064 


1.56±0.60 


0.09 


0.38±0.07 


0.03 


(1.0±0.44)xl0^° 


6.5 XlO"'^ 


< 1.50 


2 


PG0921+525 


0.035 


< 0.59 




< 0.05 




<8.6xl0*'® 


2.2x10"^ 






PG0923+201 


0.190 


< 0.35 




< 0.29 




<9.0xl0^'' 


5.9x10^" 






PG0923+129 


0.029 


9.73±2.28 


0.28 


2.42±0.13 


0.08 


(1.3±0.51)xl0^'' 


1.6x10"^^ 






PG0934+013 


0.050 


2.86±0.60 


0.26 


0.74±0.05 


0.08 


(1.2±0.48)xl0^'' 


1.8x10"^^ 






PG0946+301 


1.216 


< 0.47 




< 0.11 




<1.8xl0^^ 


1.7x10^^ 






PG0947+396 


0.205 


< 0.38 




< 0.18 




<1.2xl0^^ 


5.0x10^" 






PG0953+414 


0.234 


< 1.39 




< 0.20 




<3.8xl0^^ 


7.8x10^° 






PGlOOl+054 


0.160 


< 0.38 




0.17±0.03 


0.03 


(3.8±1.66)xl0^'' 


2.7x10 " 






PG1004+130 


0.240 


< 0.58 




0.20±0.05 


0.02 


(1.3±0.62)xl0^^ 


6.2x10^" 






PGlOll-040 


0.058 


< 0.56 




0.50±0.04 


0.03 


(l.l±0.44)xl0 


3.6x10"" 






PG1012+008 


0.186 


< 0.61 




< 0.09 




<8.2xl0^*' 


3.6x10^" 






PG1022+519 


0.044 


4.22±0.74 


0.44 


1.32±0.08 


0.19 


(1.8±0.73)xl0^'' 


1.4x10 






PG1048+342 


0.167 


< 0.33 




< 0.04 




<2.1xl0^'' 


1.3x10^" 






PG1048-090 


0.344 


< 0.33 




< 0.06 




<2.8xl0^^ 


5.4x10^" 






PG1049-005 


0.359 


1.17±0.38 


0.07 


0.17±0.07 


0.01 


(3.4±1.93)xl0^^ 


2.2x10^^ 






PGllOO+772 


0.311 


< 1.04 




0.29±0.08 


0.04 


(4.1±1.99)xl0-'--'- 


1.0x10-'-' 






PG1103-006 


0.423 


< 0.18 




< 0.09 




<3.6xl0^^ 


1.3x10-'^ 






PG1114+445 


0.143 


< 0.40 




< 0.11 




<5.2xl0"'"'-' 


4.4x10-"-' 






PG1115+407 


0.154 


2.55±0.33 


0.28 


0.46it0.03 


0.08 


(1.1±0.46) X 10-'--'- 


2.1x10-"' 






PG1116+215 


0.176 


< 3.32 




< 0.25 




<2.3x 10-'-' 


l.lxlO'^'^ 






PG1119+120 


0.050 


2.26±0.89 


0.06 


0.80±0.09 


0.03 


(1.3±0.53) X 10-"' 


5.0 x 10"" 


4.50 


1 


PG1121+422 


0.225 


< 0.35 




< 0.09 




<1.2xl0-'-' 


3.0x10^" 






PG1126-041 


0.060 


< 1.23 




1.35±0.36 


0.04 


(3.6±1.75)xl0-"' 


1.6x10-'" 


< 2.60 


1 


PG1149-110 


0.049 


< 0.64 




< 0.10 




<3.5xl0''^ 


2.2x10"" 






PG1151+117 


0.176 


< 3.30 




< 0.48 




<5.2xl0-'-' 


1.9x10^" 






PG1202+281 


0.165 


1.41±0.48 


0.14 


0.37±0.05 


0.04 


(1.0±0.44)xl0-'-' 


2.5x10-'" 


< 2.40 


1 


PG1211+143 


0.080 


< 1.82 




< 0.15 




<1.9xl0^'' 


2.8x10-'" 


< 1.50 


2 


PG1216+069 


0.331 


< 0.34 




< 0.05 




<1.9xl0^^ 


8.9x10^" 






PG1226+023 


0.158 


< 2.16 




< 0.16 




<1.0xlO-'-' 


3.7x10^^ 






PG1229+204 


0.063 


< 0.54 




0.38±0.12 


0.02 


(9.8±4.94) X lO"^ 


7.0x10 


2.40 


2 


PG1244+026 


0.048 


1.76±0.86 


0.14 


0.51±0.04 


0.04 


(7.0±2.86)xl0'-'^ 


1.8xl0"" 






PG1259+593 


0.477 


< 0.16 




< 0.04 




<3.8xl0^^ 


2.7x10-'-' 






PG1302-102 


0.278 


< 0.51 




< 0.14 




<3.7xl0^^ 


1.0x10^^ 






PG1307+085 


0.155 


< 3.47 




< 0.43 




<3.2xl0^^ 


2.6x10^" 






PG1309+355 


0.184 


< 3.17 




< 0.36 




<3.9xl0^^ 


4.4x10-'" 


< 0.61 


3 


PG1310-108 


0.034 


2.40±0.86 


0.11 


0.18±0.03 


0.01 


(1.0±0.44) X lO"'' 


1.3x10"'' 






PG1322+659 


0.168 


0.72±0.30 


0.07 


0.20±0.02 


0.03 


(5.3±2.20) X 10-"-' 


2.9x10-'" 






PG1341+258 


0.087 


0.45±0.21 


0.06 


0.11±0.02 


0.02 


(5.3±2.39)xl0'"^ 


4.7x10 






PG1351+236 


0.055 


7.54±1.05 


0.87 


2.75±0.12 


0.44 


(6.7±2.71)xl0^'^ 


1.6x10"" 






PG1351+640 


0.088 


3.12±6.54 


0.09 


1.29±0.15 


0.03 


(9.3±3.89)xl0^'' 


2.2x10^" 


4.00 


2 


PG1352+183 


0.152 


<14.14 




< 2.60 




<2.4xl0^^ 


1.7x10^" 






PG1354+213 


0.300 


< 0.27 




< 0.06 




<1.8xl0i^ 


4.2x10^" 






PG1402+261 


0.164 


< 1.59 




< 0.22 




<1.6xl0^^ 


6.8x10^" 


< 2.00 


1 


PG1404+226 


0.098 


0.88±0.37 


0.14 


0.25±0.02 


0.05 


(1.7±0.71)xl0^'' 


5.1x10 




2 


PG1411+442 


0.089 






0.31±0.04 


0.01 


(1.8±0.74)xl0^'' 




< 1.80 


2 


PG1415+451 


0.113 


1.67±0.30 


0.14 


0.86±0.06 


0.10 


(l.l±0.43)xl0^^ 


in 

1.3x10 " 


3.30 


1 


PG1416-129 


0.129 


< 0.56 




< 0.15 




<5.8xl0^*' 


8.5x10"^ 








U.OOD 






<^ u.uo 




1 V 1 fill 


1 1 V 1 nil 






PG1426+015 


0.086 


1.19±0.64 


0.03 


0.31±0.06 


0.01 


(1.7±0.73)xl0i° 


2.4x10^° 


3.60 


2 


PG1427+480 


0.221 


< 0.28 




< 0.03 




<3.4xl0i° 


2.5x101° 






PG1435-067 


0.126 


< 0.44 




< 0.19 




<4.3xl0i° 


1.7x1010 






PG1440+356 


0.079 


6.74±2.89 


0.20 


2.27±0.13 


0.10 


(1.3±0.53)xl0" 


2.0x1010 


9.00 


2 


PG1444+407 


0.267 


0.38±0.28 


0.03 


< 0.15 




(7.8±6.19)xl0i° 


8.5xlOio 


0.71 


3 


PG1448+273 


0.065 


1.98±0.59 


0.11 


0.94±0.06 


0.07 


(3.0±1.22)xl0i° 


5.7xl008 






PG1501+106 


0.036 


< 1.70 




< 0.38 




<7.9xl0°8 


4.2x1009 






PG1512+370 


0.370 


< 0.22 




< 0.07 




<3.1xl0" 


8.8x10^0 






PG1519+226 


0.137 


0.59±0.21 


0.04 


0.21±0.02 


0.02 


(3.3±1.37)xlOio 


2.7x1010 
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TABLE 1 



Continued 



source 


Redshift 


F(7.7/im) 


EW(7.7/im) 


F(11.3//m) 


EW(11.3//m) 


iSFIR 


L^—Qiim 


/I ^ 

(1) 


/ON 

(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


PG1534+580 


0.029 


1.45+0.72 


0.05 


0.44+0.08 


0.02 


(2.0+0.88) xlO"^ 


1.5x10"^ 


PG1535+547 


0.038 


0.62+0.22 


0.02 


0.08+0.03 


0.01 


(5.6+2.87) X 10"* 


3.1x10"^ 


PG1543+489 


0.399 


< 0.34 




< 0.26 




<6.3xl0^^ 


2.4x10-'--'- 


PG1545+210 


0.264 


< 1.75 




< 0.17 




<4.1xl0^i 


5.8x10^" 


PG1552+085 


0.119 


< 0.30 




0.11+0.02 


0.02 


(l.l+0.46)xl0i" 


l.OxlOi" 


PG1612+261 


0.130 


< 0.46 




0.38+0.22 


0.03 


(5.8+4.17) X 10^" 


2.1x10^" 


PG1613+658 


0.129 


3.02+1.87 


0.08 


0.77+0.09 


0.03 


(1.3+0.53) xlO^^ 


5.5x10^° 


PG1617+175 


0.112 


< 0.45 




< 0.48 




<3.2xl0^° 


1.9x10^° 


PG1626+554 


0.133 


< 0.47 




< 0.09 




<3.2xl0i° 


1.3x10^° 


PG1634+706 


1.334 


< 0.52 




< 0.11 




<2. 4x10^3 


9.4x10°^ 


PG1700+518 


0.292 


< 5.70 




< 0.20 




<6.5xl0ii 


3.2x10^1 


PG1704+608 


0.371 


< 1.04 




< 0.11 




<6.5xl0i^ 


2.6x10^1 


PG2112+059 


0.466 


< 0.24 




0.27+0.05 


0.02 


(1.2+0.52) xlO^^ 


4.8x10^1 


PG2130+099 


0.062 


4.20+1.29 


0.06 


0.55+0.21 


0.01 


(1.5+0.83) xlO^" 


2.1 xlO^" 


PG2209+184 


0.070 


1.32+0.37 


0.20 


0.29+0.03 


0.06 


(9.1+3.74) xIQOS 


3.0x1008 


PG2214+139 


0.065 


< 0.81 




< 0.28 




<1.7xl0i° 


1.6x10^0 


PG2233+134 


0.325 


< 1.44 




< 0.15 




<6.6xl0i^ 


9.4xlOio 


PG2251+113 


0.325 


< 0.55 




< 0.26 




<6.3xl0ii 


1. 5x10^1 


PG2304+042 


0.042 


< 0.46 




< 0.05 




<5.6xl0°^ 


8.8x10°* 


PG2308+098 


0.433 


< 0.27 




< 0.06 




<5. 4x10^1 


1.4x1011 


PG2349-014 


0.174 


< 0.47 




0.41+0.10 


0.05 


(1.3+0.61) x 10^^ 


4.6x101" 


2MASSJ000703. 61+155423. 8 


0.114 


3.00+0.71 


0.32 


1.00+0.10 


0.14 


(1.3+0.52) xlO^^ 


9.7x10°^ 


2MASSJ005055. 70+293328.1 


0.136 


1.76+0.34 


0.19 


0.33+0.09 


0.05 


(5.6+2.71) xlO^" 


1.5xl0i° 


2MASSJ010835. 16+214818.6 


0.285 


< 1.25 




< 0.21 




<6.6xl0^^ 


l.lxioii 


2MASSJ015721. 05+171248. 4 


0.213 


2.02+0.46 


0.33 


0.58+0.12 


0.16 


(3.5+1.56) X 10^1 


3.0x101" 


2MASSJ022150. 60+132741.0 


0.140 


< 3.31 




< 0.39 




<2.1xlO'^'^ 


2.5x101" 


2MASS,T023430. 64+243835. 5 


0.310 


< 1.16 




< 0.34 




<1.3xl0^^ 


6.4xl0i" 


2MASSJ034857.64+125547.3 


0.210 


< 1.56 




< 0.33 




<5.2xl0^^ 


2.3xl0ii 


2MASSJ091848. 63+211717.1 


0.149 


< 1.27 




0.45+0.25 


0.04 


(1.0±0.69)xl0ii 


3.2x101" 


2MASSJ095504. 56+170556.1 


0.139 


< 1.10 




< 0.28 




<1.4xlO^^ 


9.5x10"" 


2MASSJ102724. 95-121920. 4 


0.231 


< 1.49 




< 0.36 




<7.1xl0^^ 


7.0xl0i" 


2MASSJ105144. 25+353930. 7 


0.158 


< 0.99 




< 0.24 




<1. 7x10^1 


l.OxlOi" 


2MASSJ125807.46+232921.5 


0.259 


1.54+0.78 


0.09 


< 0.16 




(3.9±2.29)xl0i^ 


9.4xl0i" 


2MASSJ130005. 35+163214. 8 


0.080 


< 4.82 




< 1.05 




<1.7xl0ii 


2.2xl0i" 


2MASSJ130700. 66+233805.0 


0.275 


9.27+1.58 


0.57 


< 0.27 




(3.9+1.33)xl0i^ 


2.2xl0ii 


2MASSJ140251. 22+263117. 5 


0.187 


< 1.46 




< 0.58 




<4.4xl0^^ 


2.6xl0i" 


2MASSJ145331. 51+135358. 7 


0.139 


10.31+2.81 


0.51 


1.10+0.70 


0.16 


(2.3+1.76) X 10^1 


4.0x101" 


2MASSJ150113. 21+232908. 3 


0.258 


< 1.31 




0.19+0.08 


0.03 


(1.4+0.82)xl0^"'- 


5.0xl0i" 


2MASSJ151653. 24+190048. 4 


0.190 


< 2.22 




< 0.57 




<7.2xl0^^ 


1.6xl0ii 


2MASSJ163700. 22+222114.0 


0.211 


2.81+0.61 


0.60 


0.51+0.04 


0.15 


(2.7±1.10)xl0^^ 


2.0xl0i" 


2MASSJ165939. 77+183436. 9 


0.170 


3.43+1.21 


0.17 


0.65+0.16 


0.04 


(2.2±1.01)xl0^^ 


4.7x10"' 


2MASS,I171442. 77+260248. 5 


0.163 


1.19+0.34 


0.20 


0.32+0.07 


0.08 


(8.0+3.60) xlO^" 


1.7xl0i" 


2MASSJ222202. 22+195231. 5 


0.366 


< 0.97 




< 0.09 




<5.0xl0^^ 


2.7xl0ii 


2MASSJ222221. 12+195947. 4 


0.211 


< 1.05 




< 0.14 




<1.8xl0^^ 


4.3x10"' 


2MASSJ222554. 27+195837.0 


0.147 


1.97+0.33 


0.22 


< 0.21 




(1.2+0.40) xlO^'^ 


1.6x10"' 


2MASSJ234449. 57+122143.4 


0.199 


< 1.25 




< 0.15 




<1.8xl0^^ 


3.5x10"' 


3C6.1 


0.840 






< 0.05 




<2. 5x10^2 




3C15 


0.073 


< 0.52 




< 0.03 




<2.4xl0'-''-^ 


l.OxlO"^ 


3C20 


0.174 


< 1.20 




< 0.24 




<2.1xl0^^ 


4.2xl0°8 


3C22 


0.936 






< 0.03 




<1. 5x10^2 




3C28 


0.195 


< 0.30 




< 0.06 




<5.3xl0^'' 


8.6x10°* 


3C29 


0.045 


< 0.75 




< 0.02 




<6.6xl0°* 


2.3x10°* 


3033 


0.059 


< 0.60 




< 0.15 




<9.1xl0°^ 


2.7x10°^ 


3C33.1 


0.180 


< 1.62 




< 0.39 




<4.0xl0ii 


1.2xl0i° 


3047 


0.425 


< 0.28 




< 0.06 




<4.1xl0ii 


l.lxlQii 


3048 


0.367 


< 4.55 




< 0.53 




<4.0xl0i^ 


2.5xl0ii 


3055 


0.734 


< 0.18 




< 0.07 




<1.7xl0i^ 


8.9xl0i° 


3061.1 


0.187 


< 0.23 




< 0.05 




<4.2xl0i° 


8.3x10°* 


3065 


1.176 


< 0.42 




< 0.08 




<1. 2x10^3 


8.1x10°^ 


3075 


0.023 


< 0.33 




< 0.01 




<7.6xl00^ 


2.6x10°^ 


3076.1 


0.032 


< 0.43 




< 0.05 




<6.3xl0°^ 


9.4x10°''' 


3079 


0.255 


< 1.02 




< 0.17 




<3.9xl0" 


3.0x101° 


3083.1 


0.025 


< 0.36 




0.16+0.02 


0.09 


(4.4+1.84) xl0°8 


4.6x10°* 


3084 


0.017 






4.11+1.19 


0.02 


(7.1+3.53) xl009 


1.6xl0°9 


3O109 


0.305 


< 1.75 




< 0.26 




<1. 1x1012 


2.2x1011 


30123 


0.217 


< 0.61 




< 0.05 




<5.6xl0i° 


1.3xl0°9 


30129 


0.020 


< 0.36 




0.07+0.01 


0.06 


(1.3+0.56) x 10°* 


1.2x10°* 


30138 


0.759 






< 0.03 




<1.0xl0i2 




30147 


0.545 






< 0.05 




<6.6xl0" 




30153 


0.276 


< 0.40 




< 0.04 




<8.2xlOio 


4.7x10°* 



ScoAF 
(9) 



Ref 
(10) 



< 0.96 



8.50 



4.30 
1.60 

3.20 



2.00 
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TABLE 1 — Continued 



source 


Redshift 


F(7.7//m) 


EW(7.7/um) 


F(11.3//m) 


EW(11.3/im) 




L^ — Qiirn 


ScoAV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C172 


0.519 


< 0.17 




< 0.05 




<5.8xl0i^ 


5.5x10''^ 




3C173.1 


0.292 


< 0.31 




< 0.03 




<6. 6x101" 


1.6x10''^ 




3C175 


0.770 


< 0.17 




< 0.03 




<1.0xl0i^ 


2.8x101-*^ 




3C184 


0.994 






< 0.06 




<5.7xl0i^ 






3C192 


0.059 


< 0.43 




< 0.05 




<2.8xl0''^ 


2.1x10°* 




3C196 


0.871 


< 0.14 




< 0.04 




<2.1xl0i^ 


3.3x10^1 




3C200 


0.458 






< 0.09 




<9. 4x10^1 






3C216 


0.670 


< 0.28 




< 0.12 




<2.3xl0i^ 


2.5x10^1 




3C219 


0.174 


< 0.30 




< 0.10 




<5.9xl0i° 


6.0x10°^ 




3C220.1 


0.610 


< 0.29 




< 0.10 




<1. 8x10^2 


1.7x101° 




3C220.3 


0.680 


< 1.05 




< 0.03 




<6.9xl0ii 


1.6xl0i° 




3C234 


0.184 


< 0.90 




< 0.24 




<2. 4x10^1 


1.2x10^1 




3C244.1 


0.428 


< 0.21 




< 0.04 




<3.3xl0ii 


2.5xl0i° 




3C249.1 


0.311 


< 1.04 




0.29±0.10 


0.04 


(4.1±2.16)xl0ii 


1.0x10^1 




3C263 


0.646 


< 0.14 




< 0.07 




<9.4xl0ii 


3.6x1011 




3C263.1 


0.824 


< 0.10 




< 0.16 




<1.2xl0i^ 


1.6xl0i° 




3C265 


0.811 


0.62±0.23 


0.24 


< 0.30 




(3.4±1.62)xl0i^ 


2.6xl0ii 




3C268.1 


0.970 


< 0.15 




< 0.08 




<3.2xl0i2 


2.2xl0i° 




3C270 


0.007 






0.60±0.04 


0.09 


(1.4±0.56)xl0°^ 


5.4x10°''' 




3C272 


0.944 






< 0.02 




<1. 1x10^2 






3C272.1 


0.003 






1.70±0.12 


0.33 


(2.0±0.04)xl0°^ 


3.0x10°''^ 




3C273 


0.158 


< 2.16 




< 0.16 




<1. 0x10^1 


3.7xl0ii 




3C274 


0.004 






< 0.97 




<2.3xl0°8 


4.3x10°'' 


< 11.7 


3C274.1 


0.422 


< 0.19 




< 0.08 




<3. 7x10^1 


3.0x10°^ 




3C275.1 


0.555 


< 0.15 




0.09±0.02 


0.09 


(5.1±2.29)xl0ii 


5.2x101'' 




3C280 


0.996 


< 0.09 




< 0.09 




<1.9xl0i^ 


2.3x1011 




3C292 


0.710 






< 0.09 




<3.4xl0i^ 






3C293 


0.045 


3.96±0.70 


0.62 


1.27±0.10 


0.41 


(1.7±0.71)xl0i'' 


9.2x10"*^ 




3C295 


0.464 


< 0.13 




< 0.24 




<3. 1x10^1 


3.5x10°-' 




3C298 


1.436 


< 0.30 




< 0.07 




<1.7xl0i^ 


1.2xl0i^ 




3C300 


0.270 


< 0.34 




< 0.06 




<1. 2x10^1 


1.5x10°" 




3C303.1 


0.267 


< 0.38 




0.09±0.02 


0.18 


(6.9±3.21)xl0"' 


5.3x10°'-' 




3C309.1 


0.905 


< 0.11 




< 0.03 




<1.9xl0i^ 


3.3xl0ii 




3C310 


0.053 


< 0.30 




< 0.03 




<1.3xl0'''5 


1.5x10°* 




3C315 


0.108 


< 0.44 




0.17±0.02 


0.50 


(1.4±0.61)xl0i'' 


5.4x10°* 




3C318 


1.574 


< 0.51 




< 0.07 




<2.3xl0i^ 


2.8xl0ii 




3C319 


0.192 


< 0.23 




< 0.08 




<5.6xl0"' 






3C321 


0.096 


6.51±1.04 


0.49 


< 0.28 




(1.7±0.57)xl0ii 


6.1x10°" 


< 4.70 


3C323.1 


0.264 


< 1.75 




< 0.17 




<4.1xl0ii 


5.8xl0i° 




3C325 


1.135 


< 0.10 




< 0.04 




<3.3xl0i^ 


9.3xl0i° 




3C326 


0.089 


< 0.62 




< 0.11 




<1.7xl0i'' 


3.2x10°* 




3C330 


0.550 


0.25±0.07 


0.29 


< 0.02 




(3.8±1.54)xl0ii 


2.8xl0i° 




3C334 


0.555 


0.58±0.21 


0.17 


< 0.03 




(l.l±0.50)xl0i^ 


1.5xlOii 




3C336 


0.927 






< 0.08 




<6.6xl0i^ 






3C337 


0.635 






< 0.05 




<1.2xl0i^ 






3C340 


0.775 






< 0.03 




<l.lxl0i^ 






3C343 


0.988 






< 0.04 




<3. 2x10^2 






3C343.1 


0.750 






< 0.02 




<7.2xl0i'^ 






3C348 


0.154 


< 0.81 




< 0.19 




<1. 1x10^1 


8.3x10°* 




3C351 


0.371 


< 1.04 




< 0.11 




<6. 5x10^1 


2.6xl0ii 




3C352 


0.806 






< 0.05 




<2.3xl0i^ 






3C356 


1.079 


< 0.17 




< 0.06 




<5.2xl0i^ 


8.0x101° 




3C371 


0.051 


< 2.12 




< 0.13 




<5.1xl0°^ 


8.6x10°" 




3C380 


0.692 


< 0.17 




< 0.09 




<1.3xl0i^ 


3.5xl0ii 




3C381 


0.160 


< 0.48 




< 0.05 




<3.0xl0i° 


1.5xl0i° 




3C382 


0.057 


< 0.97 




< 0.12 




<6.3xl0°^ 


1.7xl0i° 




3C386 


0.016 






< 0.04 




<1.4xl0''^ 


9.0x10°''' 




3C388 


0.091 


< 0.36 




< 0.06 




<8.8xl0°^ 


5.9x10°* 




3C390.3 


0.056 


< 0.63 




< 0.16 




<7.9xl0°^ 


8.3x10°" 


< 10.3 


3C401 


0.201 


< 0.27 




< 0.03 




<2.9xl0i° 


1.1x10°" 




3C403.1 


0.055 


< 0.25 




< 0.02 




<7.1xl0°^ 


1.7x10°* 




3C405 


0.056 


< 3.28 




< 0.55 




<3.3xl0i° 


3.5x10°" 


< 1.90 


3C427.1 


0.572 


< 0.22 




< 0.03 




<5.6xl0" 


3.4x10°" 




3C433 


0.101 


< 0.83 




< 0.22 




<4.9xl0i° 


2.0x101° 




3C436 


0.214 


< 0.40 




< 0.05 




<5.4xlOio 


1.1x10°" 




3C438 


0.290 


< 0.35 




< 0.04 




<9.3xl0i° 


1.5x10°" 




3C441 


0.708 






< 0.06 




<1.9xl0i2 






3C445 


0.056 


< 1.48 




< 0.30 




<1.7xl0i° 


1.7x101° 




3C452 


0.081 


< 0.50 




< 0.07 




<7.9xl0''^ 


2.7x10°" 




3C465 


0.030 


< 0.86 




< 0.21 




<2.6xl009 


3.5x10°* 





Ref 
(10) 
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TABLE 1 — Continued 



source 
(1) 


Rcdshift 
(2) 


F(7.7Mm) 
(3) 


EW(7.7Atm) 
(4) 


F(11.3/im) 

(5) 


EW{11.3Atm) 

(6) 


(7) 


Ls-e^m Sco^^ 
(8) (9) 


Ref 
(10) 





Note. — Column (1): Sources. Column (2): Rcdshift. Column (3): The observed-frame 7.7 fim aromatic flux in the unit of 10 erg s ^ 
cm~^. Column (4): The rest-frame EW of 7.7 fim PAH in the unit of /j.m. Column (5): The observed-frame 11.3 (J,m aromatic flux in the unit of 
10~^^ erg s"^ cm~^. Column (6): The rest-frame EW of 11.3 fim PAH in the unit of fim. Column (7): The star-forming IR luminosity in the 
unit of Lq. Column (8): The mid-IR luminosity in the unit of Lq integrated from 5 to 6 fim. A factor of 22.6 can be applied to convert it to the 
total IR luminosity (3-1000 fim) based on the quasar template of lElvis etlill lll994r) . Column (9): The CO flux in the unit of Jy km s ^. Column 

(10): Reference for column (9) ■ 

REFERENCES: (l lEvans et al.l | |2001|) : (2 lScoville et al.l II2003D : O lCasoU fc LoinardI | |2001|) : (4 1Scoville et al.l J1993D : (S lEvans et al.l | |2005D 



TABLE 2 

The Star Formation Fraction at Three MIPS Bands as a Function of 
THE MiD-IR Luminosity 



MIPS band 


Q 


13 


Correlation 


A1I(MIPS 24 iim) 
A11(MIPS 70 urn) 
A11(MIPS 160 Aim) 


0.6± 1.3 
1.2± 1.0 
2.3± 5.5 


-0.18±0.13 
-0.17±0.10 
-0.27±0.51 


-0.22± 0.15 
-0.32± 0.17 
-0.15± 0.29 


PC(MIPS 24 Atm) 
PG(MIPS 70 Atm) 


0.2± 1.3 
2.4± 1.6 


-0.15±0.13 
-0.29±0.15 


-0.22± 0.19 
-0.43± 0.19 


2MASS(MIPS 24 Aim) 
2MASS{MIPS 70 Atm) 


0.1±12.3 
1.2± 6.9 


-0.11±1.15 
-0.14±0.65 


O.Olit 0.43 
-0.05± 0.48 



Note. — Log{Fva.cfp^^) = a + f3xLog{LMiR) 



TABLE 3 

The Star Formation Fraction at Three MIPS Bands as a Function of 

THE EdDINGTON RATIO 



MIPS band 


a 




Correlation 


PG(MIPS 24 Aim) 


-1.6± 0.3 


-0.10±0.12 


-0.18± 0.21 


PG(MIPS 70 Aim) 


-1.3± 0.3 


-0.32±0.12 


-0.60± 0.17 



Note. — Log(Frac^i?^) = a + l3xLog(LMiR/LEdd) 



TABLE 4 

Fractional Bivariate Luminosity Function for PC quasars 







Mb (mag) 




Log(ifi,^[Lo]) 


-25.83 


-24.57 


-23.31 


-22.05 


10.06 
10.75 
11.43 
12.11 


0/0± 1.00 
0/0± 1.00 
l/3± 0.38 
1/14± 0.07 


0/0± 1.00 
2/2± 1.00 
2/8± 0.20 
1/17± 0.06 


3/6± 0.35 
7/19± 0.16 
3/26± 0.07 
0/26± 0.00 


8/12± 0.30 
5/14± 0.19 
1/15± 0.07 
0/15± 0.00 



TABLE 5 

Fractional Bivariate Luminosity Function for 2MASS quasars 



^og{L[J/j,[LQ]) 




Log{LK [L©]) 




10.88 


11.26 


11.64 


11.06 


4/5± 0.54 


3/2± 1.37 


0/0± 1.00 


11.70 


l/9± 0.12 


2/8± 0.20 


0/0± 1.00 


12.33 


0/10± 0.00 


1/12± 0.09 


0/2± 0.00 



24 



Shi et al. 



TABLE 6 

Fractional Bivariate Luminosity Function for SCR radio galaxies 

AND QUASARS 





Log(Li5iMHz[W Hz 


■isr-i]) 


Log(Lfi/^[Lo]) 


25.13 


26.23 


27.32 


10.54 
12.01 


4/9± 0.27 
0/10± 0.00 


l/3± 0.38 
1/21± 0.05 


0/l± 0.00 
0/9± 0.00 



TABLE 7 

Best-fitting parameters to star-forming IR LF of PG quasars 



Object 


Log((^*[Mpc ^ mag ^]) 


Log{L*[Lc3]) 


a 


PG(Mb <-21) 


-7.88± 0.29 


11.45± 0.17 


-1.18± 0.24 


PG(Mb <-23) 


-8.37± 0.35 


11.49± 0.42 


-0.28± 1.49 



Note. — The formula of luminosity function is a Schechter function: ^{L)dL = (-^)"exp(— -^)-^. 



